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Abstract 

When small RNAs are loaded onto Argonaute proteins they can form the 
RNA-induced silencing complexes (RISCs), which mediate RNA interference. 
RISC-formation is dependent on a shared pool of Argonaute proteins and 
RISC loading factors, and is thus susceptible to competition among small 
RNAs for loading. We present a mathematical model that aims to under- 
stand how small RNA competition for the PTR resources affects target gene 
repression. We discuss that small RNA activity is limited by RISC-formation, 
RISC-degradation and the availability of Argonautes. Together, these obser- 
vations explain a number of PTR saturation effects encountered experimen- 
tally. We show that different competition conditions for RISC-loading result 
in different signatures of PTR activity determined also by the amount of 
RISC- recycling taking place. In particular, we find that the small RNAs 
less efficient at RISC-formation, using fewer resources of the PTR pathway, 
can perform in the low RISC-recycling range equally well as their more ef- 
fective counterparts. Additionally, we predict a novel signature of PTR in 
target expression levels. Under conditions of low RISC- loading efficiency and 
high RISC-recycling, the variation in target levels increases linearly with the 
target transcription rate. Furthermore, we show that RISC-recycling deter- 
mines the effect that Argonaute scarcity conditions have on target expression 
variation. Our observations taken together offer a framework of predictions 
which can be used in order to infer from experimental data the particular 
characteristics of underlying PTR activity. 
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Introduction 

Post-transcriptional regulation (PTR) is mediated by RNA-induced si- 
lencing complexes (RISCs) assembled from Argonaute proteins and microRNA 
(miRNA) or small- interfering RNA (siRNA) molecules^. RISCs act on mRNA 
transcripts via the RNA-interference pathway increasing the mRNA turnover 
(degradation) rate via cleaving^ or destabilization^. A single RISC can act 
multiple times like an enzyme^, even when target destabilization rather than 
cleaving takes place mRNA destabilization and cleaving differ in the 
number of times a RISC can operate before dissociation or degradation of 
the loaded small RNA. For example, let-7 miRNA-programmed RISCs were 
found to catalyze each approximately 10 target molecules^. Another study 
showed that let-7 siRNA-programmed RISCs participated each in 50 rounds 
of catalysis on average^. However, experiments using miR223-programmed 
RISCs constitutively expressed in human 293T cells reported an average of 
only two rounds of activity per RISC'^. Conceptually, it is simpler to study 
PTR in the two limits of RISC operational activity^'^'^"'^^: the non-enzymatic 
stoichiometric mode, where each RISC degrades with the transcript it has 
targeted, and the fully enzymatic catalytic mode, where RISC-degradation 
is independent of mRNA-targeting. Part of the scope of this work is to 
understand the dynamics of PTR activity in these two distinct operational 
regimes. 

Post-transcriptional regulation is a mediated process and as such, is sub- 
ject to competition and saturation effects. Mature small RNAs need to be 
loaded onto Argonautes in order to be activated. This is a process which in- 
volves competition to access and be loaded to the Argonautes and can lead to 
saturation conditions for the underlying RISC-forming machinery of the cell. 
This saturation can take place either through the large number of mature 
small RNAs present, or for example under Argonaute scarcity conditions. 

The experimental evidence for competition and saturation effects in PTR 
comes from multiple sources. For example, it was shown that highly stable 
chemically modified siRNAs are able to outcompete the less stable native 
siRNAs^^. Another study presented in vitro evidence of intra-siRNA com- 
petition and also siRNA-miRNA competition for Argonautes Addition- 
ally, up-regulation of target gene expression levels observed in transfection 
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experiments was shown to be likely due to RISC saturation effects Reg- 
ulation during the maturation process of small RNAs^^, miRNA sponges 
and RISC co- factors can also induce competition-like and saturation-like 
effects of various strengths to different small RNAs. On the other end of the 
competition spectrum, with Argonautes highly abundant, it has been found 
that small RNAs compete to access targets and that transcript abundance 
dilutes the small RNA activity This is because transcripts with target sites 
for the same miRNA can compete for binding, in which case the titration of 
miRNAs by a certain target can release the repression from other targets 

In the first part of this paper we construct a mathematical model of 
miRNA biogenesis and activity which explicitly takes into account RISC- 
formation. Understanding how saturation and competition in RISC- formation 
can influence the dynamics of PTR is the second aspect of this work. We 
will show that there is a distinctively different signature in the variation 
of target transcript expression depending on whether the mode of PTR is 
stoichiometric or catalytic. 



Methods 

Synthesizing a minimal model of PTR. We construct a model of 
post transcriptional regulation based on the loading of small RNAs onto 
Argonautes. Our terminology is equally applicable to miRNA and siRNAs, 
and comparisons with experiments will draw on both prokaryotic and eu- 
karyotic model systems. The biogenesis of a single species of small RNA, the 
expression of targets, and RISC activity is modeled along the lines of pre- 
vious work®'^'^. Creation and degradation of molecules of a particular small 
RNA species (m), the target mRNA (R), and the target protein (P) are 
described by the following set of reactions 

^ m, (small RNA) (1) 

dill 

^ R, (target mRNA) (2) 



R%R+P 
P A 



(target Protein). (3) 



The values for all parameters are provided in Table 1. We denote the number 
of small RNA bound to Argonautes by Am-, and the number of unbound 
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Argonautes by A. For simplicity we take the total number of Argonautes 
Aq = A + Ara to be constant. We then describe RISC assembly at rate /c^ 
and RISC disassociation at rate (Ia by 

m + A ^ (RISC). (4) 

We stress that the RISC- format ion process is independent of the RISC tar- 
geting activity, the latter described by the reaction 

A^ + r'^^\^ stoichiometric 1 

1^ Am catalytic J ^ ^ ^ ^ 

with a rate djim., and an outcome which depends on the mode of PTR op- 
eration, stoichiometric or catalytic. For simplicity, we describe single-seed 
targets only, but our result are not conditional on this assumption. Accord- 
ing to Eq. (5), a targeted mRNA via either mode of PTR cannot be released 
back into the pool of active mRNAs available for translation. We consider 
the set of Eqs. (1-5) as a minimal model of PTR mediated by small RNAs. 

Under this minimal model, the number of Argonautes available is deter- 
mined by the parameter Aq. This parameter is affected by Argonaute biogen- 
esis and turnover, as well as other small RNA species binding to Argonautes, 
so it incorporates an element of the competition faced by the modeled small 
RNAs to access the Argonautes. The actual formation of RISCs is determined 
by the rate kA, which depends primarily on the individual efficiency that the 
small RNA species has to load to an Argonaute. The abundance of avail- 
able RISC loading factors, which in turn depends on the presence of other 
competing small RNA species, can influence /c^ as well. Overall, the ratio 
Mq = dA/kA is a good measure of the effectiveness of small RNAs at RISC- 
formation which also takes into consideration how stable the small RNAs are 
when incorporated in RISCs. Low values of Mq indicate that the small RNAs 
considered are highly effective at incorporating and stabilizing themselves in 
RISCs^. High values of Mq reduce the inherent RISC-incorporating effec- 
tiveness, or increase the RISC-instability of the modeled small RNAs^^ The 
presence of other equally potent or highly expressed small RNAs^^, or the 
saturation of co- factors involved in RISC-assembly^^, are additional elements 
that contribute to high Mq values. 

In our study, the Fano factor of the number of molecules is the measure 
used to capture the strength of variation in different molecular species. It is 
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defined as the ratio of tlie variance over tlie mean of the expression level of 
a molecular species. For example, a Fano factor of target transcripts equal 
to unity arises from mRNA expression and degradation at constant rates 
A Fano factor exceeding unity indicates additional sources of mRNA varia- 
tion^^. By harvesting repeated measurements under steady-state conditions, 
we collect sufficient statistics of the PTR module simulated using the Gille- 
spie algorithm (direct method) In order to draw direct comparisons with 
previous theoretical work^°, the regime of PTR operation we present here 
is that of Salmonella. However, PTR in mammals was simulated as well, 
verifying that our results are not limited to PTR in bacteria only but extend 
over a broader scope. 

Results 

RISC-formation introduces an upper bound to the efficiency of 
small RNAs. The system of reactions described by Eqs. (1-5) is approxi- 
mated in the limit of large copy-numbers^^ by the following rate equations 



dm 



dt 



km - dm.m - kAmA + dA^m, (6) 



dt ~dt ^ kAmA - dAA^ -durn-AmR : (7) 

stoichiometric 

kR — duR — durnAmR, (8) 



dR 

dp 

~dt 



kpR-dpP. (9) 



The term over the bracket is present only in stoichiometric regulation. If 
correlations among the constituents are neglected, an analytic solution for 
Eqs. (6-9) can be found at steady-state. The part of the solution associated 
with RISC-formation, when combined with our assumption of A + A^ = Aq, 
yields 

m/Mo , , 

stoichiometric 
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For the constituents involved in RNA interference we find 




stoichiometric 



(11) 



R 



(12) 



1 + 




where is given by Eq. (10) and in stoichiometric PTR can depend also on 
the number of mRNAs. In Eq. (10), = dA/dnm is a threshold that emerges 
in the stoichiometric regime only and controls the targeting efficiency of the 
RISC complex before it is destabilized with the target. 

We find a Michaelis-Menten functional form for Eq. (10), which is not 
surprising since the Argonautes operate in a manner similar to enzymes^. 
Experimental evidence for Eqs. (10,12) comes from mammalian cell lines. 
Cuccato et al. fit different models of PTR to data from cell lines express- 
ing EGFP or tTA proteins. Some of these models had a basis in explicit 
biochemical reactions, but the best fit was provided by a phenomenological 
model without a simple biochemical basis. These best-fitting equations were 
of the form of Eqs. (10,12) and we conclude that RISC-formation is the key 
mechanism behind their success. 

When the expression level of a particular small RNA increases in the cell, 
this species is more likely to be loaded unto Argonautes, be incorporated 
into RISCs and repress its targets. In this case, the rate-limiting factor 
in the strength of target repression is set by Mq, i.e., the RISC-assembly 
mechanism. When the small RNA expression reaches m ^ Mq, the target 
repression rate begins to level off as saturates due to competition between 
small RNAs and saturation of the RISC-assembly machinery. This effect is 
the basis behind particular experimental observations. For example, it was 
found that differences of orders of magnitude among small RNAs and target 
transcripts do not suffice to fully suppress target expression^^'^^, or that only 
moderate changes in transcriptome expression levels are observed when small 
RNAs are transfected into cells^'^'^'^'^'^'^. It is not only the amount of small 
RNAs present in the cell relative to the number of targets that determines the 
strength of transcript suppression, but also the relative competition that the 
small RNAs face in order to perform their tasks, and the number of targets 
that each small RNA has^^. As the extra term in Eq. (10) reveals, which is 
present only in stoichiometric PTR, the mode of PTR operation plays also 



6 



a crucial role in this context. In stoichiometric PTR, RISCs are degraded 
with the targets, but in fully catalytic PTR this does not occur. As a result, 
target derepression from regulation is easier under low RISC- recycling than 
high RISC- recycling conditions. 

RISC-assembly influences the variance of target transcript ex- 
pression. Different species of small RNAs differ in their efficiency of RISC- 
incorporation. In order to simplify the discussion, we present simulation 
data from two functional regimes of PTR operation: the high efficiency (HE) 
regime characterized by Mq <C m values, and the low efficiency (LE) regime 
characterized by Mq S> m values. We stress that the same small RNA species 
transfected or natively expressed in different cell types can find itself facing 
HE or LE conditions, depending on the presence and expression levels of 
other small RNAs, and on the amount of freely available Argonaute proteins. 
Furthermore, the mode of PTR can be affected by competition as well. For 
example, it was found in mammalian cells that miRNAs are randomly sorted 
to the Argonautes, although the only Argonaute in mammals with RNA sheer 
activity is Ago2^^. As a result, the same small RNA when loaded onto an 
Ago2-programmed RISC will perform cleaving-mediated PTR activity faster 
and possibly with higher RISC- recycling compared to Agol,3,4-programmed 
RISCs, which rely on the action of nucleases and have in principle lower 
RISC- recycling. The sorting of small RNAs to the Argonautes is random^^, 
therefore it is also highly susceptible to competition conditions induced by 
the expression of other small RNAs. 

Highly efficient small RNAs saturate the PTR pathway. First, 
we simulate cell transfections with highly efficient small RNAs and under 
conditions of Argonaute abundance. Figure 1 shows the Fano factors and 
average molecular numbers per cell as functions of the target gene transcrip- 
tion rate /cr. The case of HE stoichiometric PTR is shown in Figure lA. In 
stoichiometric PTR, RISCs are degraded with the targets. However, when 
the target transcription rate is lower than the small RNAs transcription rate 
{kji < km), then RISC-target degradation has a weak effect in RISC and 
substrate copy-numbers. The highly efficient small RNAs titrate the Arg- 
onaute supply of the cell, as the inset of Figure lA shows, and additionally a 
number of free small RNAs remains as substrate. On the other hand, targets 
with transcription rates close to the small RNA production rate {kn ~ km) 
are expressed at compensatory levels to those of the small RNAs, and begin 
to alleviate the PTR-mediated repression This occurs via substrate de- 
pletion, since the newly available Argonautes from the degraded RISCs are 
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quickly recruited by the highly efficient small RNA substrate. The result is a 
very sharp transition from target repression to target derepression for those 
mRNAs with transcription rates > km- The main graph of Figure lA 
shows this sharp signature in expression level variation, especially for the 
substrate and RlSCs. Highly efficient small RNAs force RISC-target degra- 
dation to take place over a shorter k^ range, and repression drops sharply 
for targets with kf{ > k^. This property can be desirable for example for 
miRNAs or siRNAs involved in sharp activation and deactivation transitions 
of genes during the developmental phase of organisms. 

Figure IB corresponds to HE catalytic PTR, where the range of the 
target transcription rate is increased 6-fold. Here, there is no depletion of 
the free small RNA substrate taking place, and the RISC-target interaction 
does not degrade RISCs. Titration of the Argonautes by the highly efficient 
small RNAs takes place again, as the inset of Figure IB shows, but this 
time over the full range of mRNA transcription rates. The PTR module 
operates at maximum target repression capacity, and targets are strongly 
repressed and show small expression variation for any range of transcription 
rates. Overall, we notice that targets with transcription rates far from the 
derepression threshold, whether catalytically or stoichiometrically regulated, 
show variation of unity Fano factor strength. 

Similar profiles in target expression variation were reported by a theo- 
retical study which did not include the RISC-assembly process explicitly 
In view of our results and analysis so far, this is not surprising. HE PTR 
conditions are conditions of fast RISC assembly, in which case an effective 
model which does not directly include this process is accurate. However, 
it should be emphasized that HE conditions for a given small RNA gene, 
are LE conditions for the rest of the small RNA genes expressed in the cell. 
If Argonautes can be as efficiently titrated by the product of a single or a 
few highly efficient small RNA genes as Figure 1 shows, then the rest of the 
expressed small RNAs are facing Argonaute scarcity conditions. 

High RISC-recycling for less efficient small RNAs induces a 
novel signature in target expression variation. We now investigate 
post-transcriptional regulation mediated by small RNAs which have a low 
efficiency at RISC-formation. Figure 2A shows the case of stoichiometric LE 
PTR. As the free small RNA substrate is less responsive to RISC-formation, 
Argonaute titration becomes less efficient and fewer RISCs are formed, the 
latter becoming easier to be degraded with targets. The transition from tar- 
get repression to target derepression becomes smoother, and the peak Fano 
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factor value in target levels is reduced compared to the case of highly ef- 
ficient small RNAs. It should be mentioned that d^m, the rate at which 
RISCs target transcripts, is the same for both HE and LE conditions. The 
term dRm-^mR in Eq. (8) implies that target variation is affected by the RISC 
copy-number. Fast RISC-formation leads to a larger pool of RISCs. Conse- 
quently, during RISC-depletion close to the derepression threshold, targets 
experience higher expression variation. This is the reason why the peak of the 
Fano factor of target expression is lower for LE than HE small RNAs. Over- 
all, we find that in stoichiometric PTR the less efficient small RNAs, which 
are not demanding on PTR resources, can regulate targets as effectively as 
highly efficient small RNAs which tend to saturate the PTR network. Viewed 
differently, our result suggests that stoichiometric PTR involving small RNAs 
of low efficiency has a capacity for many different independent "channels" , 
where the different types of small RNAs are concurrently involved in inde- 
pendent stoichiometric PTR tasks. 

Figure 2B shows the case of catalytic LE PTR with the kn range increased 
6-fold. As expected for PTR not operating at capacity, target derepression 
occurs at lower kn values compared to HE PTR. However, a qualitatively 
different profile of PTR activity emerges. The Fano factor, and hence the 
strength of variation in target mRNA numbers increases linearly with the 
target transcription rate. Target transcripts with higher transcription rates 
than their regulators experience higher variation in their expression levels 
the higher k^ becomes. Target repression is again more effective compared 
to stoichiometric PTR, but now the profiles in variation of target expression 
under stoichiometric and catalytic LE PTR are distinctively different. 

It is known that the expression variation measured by the Fano factor for 
the biogenesis of mRNA under constant transcription and degradation rates 
is independent of the transcription rate^°. This scenario is recovered far from 
the derepression threshold k^^ k^ in stoichiometric LE PTR, and in both 
stoichiometric and catalytic HE PTR. However, we find under catalytic LE 
PTR that the Fano factor becomes linearly dependent on the target tran- 
scription rate. The root cause of this effect is suggested by Eq. (12) showing 
that target abundance is affected by RISC abundance. Due to the decreased 
ability at RISC-formation of LE small RNAs, the average copy-number of 
targets responds to sustained changes in RISC numbers. As a result, there 
are long time intervals where target turnover readjusts to changes in the 
number of RISCs. Any increase in k^ enlarges the pool of target transcripts 
-R, but also increases the turnover rate associated with the term dRrn-^mR 
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of Eq. (8), simply because the same amount of RISC is now exposed to more 
targets. Consequently, a linear increase of target expression variation with fc^ 
follows. On the other hand, small RNAs efficient at RISC-formation induce 
fast changes to RISC copy-numbers, averaging over time to a constant rate of 
RNA interference activity. This is the reason why the LE range of PTR oper- 
ation could not be captured by previous models lacking the RISC-formation 
process 

RISC-recycling determines the response of PTR to Argonaute 
scarcity conditions. We have discussed that Argonaute abundance condi- 
tions for small RNAs highly efficient in RISC-formation can result in Arg- 
onaute scarcity conditions for other less efficient but concurrently expressed 
small RNAs. Additionally, highly efficient small RNAs due to competition 
from their isoforms, or from different small RNAs of similar efficiency, can 
face also Argonaute scarcity conditions. So far, we have investigated cells 
with Aq = 500 Argonautes available. In what follows, we are going to probe 
lower values of Aq for the LE operational regime of PTR. 

Figure 3A considers stoichiometric LE PTR close to the derepression 
threshold = km- As Argonaute availability increases, more small RNAs 
can be incorporated into RISCs, but on the other hand more RISCs can 
be degraded with targets. Overall, the derepression threshold at kn = km 
remains fixed, but the average copy-number of substrate small RNAs and 
targets decreases with increasing Aq, as Eqs. (11-12) predict as well. Conse- 
quently, as Aq levels increase in the cell, m and Am show greater expression 
level variation which affects also the expression level variation of targets. 

The picture reverses for catalytic LE PTR, as shown in Figure 3B. RISCs 
can no longer be degraded with targets in catalytic PTR, which renders the 
number of RISCs independent of the number of targets. On the other hand, 
the RISC copy-number is dependent on the Argonaute copy-number, but the 
amount of variation of the RISC copy-number is independent of Aq as shown 
in Figure 3B. Furthermore, for the substrate small RNAs at steady-state, 
both the copy-number as Eq. (11) predicts, and strength of variation become 
independent of Aq. Lower Argonaute numbers reduce the PTR- mediated 
degradation rate of targets close to the basal degradation rate. Consequently, 
the lower Aq becomes, the more that variation in RISC copy-numbers influ- 
ences variation in the turnover of targets. It should be mentioned that the 
monotonic increase of the Fano factors with decreasing Aq shown in Figure 3B 
is determined by the relative strength of the PTR-mediated degradation rate 
of the target mRNA compared to the basal degradation rate. At extreme 
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Argonaute scarcity conditions, where dfimAmR ^ c?_r, a non-monotonic pro- 
file emerges, and an inflection point appears close to the crossover range 
dfim.AmR ~ dfi (data not shown). In Figure 3B, the crossover range corre- 
sponds essentially to Aq = 0, therefore no inflection point is observed. 

Under LE PTR, the impact of Argonaute scarcity conditions on target 
expression variation is dependent on RISC-recycling. In the limit of weak 
RISC-recycling, a reduction in Argonaute availability reduces the strength 
of variation in target expression, but in the limit of strong RISC- recycling 
the opposite holds true. This finding complements our conclusion from our 
study of stoichiometric regulation that less efficient small RNAs can be fa- 
vored because they require less Argonautes to perform tasks equally well. If 
multiple stoichiometric PTR "channels" are open, the channel with the least 
efficient or lowest expressed small RNAs will be operating under Argonaute 
scarcity conditions. However, as Figure 3A shows, the level of "noise" in 
PTR for this channel will be reduced compared to the rest of the activated 
channels. On the other hand, for catalytic regulation when multiple types 
of small RNAs are coexpressed, Argonaute scarcity conditions induce an in- 
crease to the strength of variation in the expression of targets, leading to 
more "noisy" PTR activity. 

Discussion 

Understanding competition effects within the PTR pathway is of central 
importance to PTR. For example, transfection protocols using synthesized 
small RNAs to target particular genes, require taking into account and opti- 
mizing for the strength of perturbation induced by the small RNA introduced 
to the cellular environment. Highly expressed small RNAs, or small RNAs 
efficient in RISC-formation face low competition for RISC-loading in the 
transfected cell lines. However, the natively expressed small RNAs in these 
cell lines face high competition, which may not be part of the physiological 
range of activity. Endogenous aberrant expression of miRNAs in cancer cells 
can lead also to PTR competition conditions'^^. Here, abnormalities in the 
expression or efficiency of small RNAs due to mutations or DNA damage can 
perturb the rest of the normally functioning PTR pathway. 

Accumulated insight into the mechanisms of PTR suggests that compe- 
tition effects are ubiquitous also under normal physiological conditions. For 
example, most mammalian miRNA genes have multiple isoforms^. When 
different isoforms have similar efficiency at RISC-formation, one can expect 
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competition to set in even when low numbers of miRNA or siRNA genes 
are expressed. Similar conditions also arise for PTR in prokaryotes, with 
differences to eukaryotic PTR being quantitative rather than qualitative'^^. 

It is important to understand that within the same cellular environment, 
in principle different PTR conditions can apply for the different small RNA 
species and their targets. If multiple highly efficient or abundant small RNA 
genes are natively expressed, their targets will exhibit the variation predicted 
under LE instead of HE PTR conditions, since the different species of small 
RNAs will compete for the Argonautes. Additionally, targets of the rest of 
the natively expressed small RNAs in these cell lines might experience PTR 
under Argonaute scarcity conditions, if the majority of the Argonautes are 
titrated by the more efficient small RNAs. 

We have developed a mathematical model of RISC-formation from which 
several principles of post-transcriptional regulation emerge. Argonaute pro- 
teins operate on the substrate of free small RNAs as enzymes, and when 
modeled as such, the functional form of Eqs. (10,12) emerges. This func- 
tional form was used as a heuristic without biochemical basis to model PTR 
data^^, where it gave a superior fit to the data compared with several other 
models not incorporating RISC-formation. Our results indicate that the up- 
per bound in the efficiency of small RNAs is set not only by the number 
of Argonautes in the cell but also by the RISC-degradation rate (Ia present 
in the RISC-saturation thresholds Mq and Rq in Eq. (10). Furthermore, 
RISC turnover is affected in stoichiometric PTR by the joint RISC-target 
degradation process and becomes dependent on as well. 

Small RNAs play a crucial role during the developmental phase of or- 
ganisms by activating or deactivating networks of genes, or by maintaining 
expression thresholds. Our mathematical model predicts that in stoichiomet- 
ric PTR the target derepression transition is sharp under HE conditions but 
smoother and less "noisy" under LE conditions. On the one hand, highly 
efficient small RNAs can perform their tasks better and faster. On the 
other hand, they can outcompete the less efficient small RNAs concurrently 
involved in PTR, thus can saturate the capacity of the PTR pathway for 
multi-tasking. The profiles of stoichiometric PTR under LE and HE condi- 
tions are similar, but small RNAs less efficient at RISC-formation use fewer 
of shared cellular resources. This opens up the possibility for greater num- 
bers of independent PTR pathways to coexist, allowing for more complex 
cellular functionality to emerge. 

For the case of small RNAs with low RISC-formation efficiency we found 
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a novel signature in the variation of target expression levels. Although tar- 
get suppression is less severe compared to PTR mediated by highly efficient 
small RNAs, our model predicts that the variation in the target expression 
levels (measured by the Fano factor) increases linearly with the target tran- 
scription rate. Transcripts targeted catalytically and kept at a threshold level 
show larger expression variation around the threshold when they are targeted 
under low efficiency than high efficiency PTR conditions. This finding sug- 
gests that for example an organism with two post-transcriptionally regulated 
paralogous genes, one with accumulated mutations in the seed-region, and 
the other perfectly complementary to the small RNA regulators, will show 
different RNA interference dynamics, depending on which copy is expressed. 
The transcripts with imperfect complementarity to the small RNA regulators 
tend to follow the stoichiometric mode of PTR and will be stably repressed 
or kept at a threshold. The transcripts with perfect complementarity follow 
the catalytic mode and when kept at a threshold will show greater expres- 
sion level variance. Essentially, this organism will possess the capacity to 
control the expression profile of a particular protein simply by activating the 
corresponding copy of the paralogous gene. 

We have presented data of PTR corresponding to two limiting RISC- 
recycling rates: stoichiometric and catalytic. However, our investigation was 
extended also to intermediate RISC-recycling rates. Highly efficient small 
RNAs of intermediate RISC-recycling were studied before For low effi- 
ciency small RNAs we find that the fully catalytic profile of Figure 2B is 
modified and an infiection point appears at a level of target expression where 
the joint RISC-target degradation is able to deplete the small RNA substrate 
(data not shown). The higher the RISC- recycling rate, the further away in 
kfi ^ km values this infiection point appears. Overall, our results and conclu- 
sions based on stoichiometric and catalytic regulation are not significantly 
modified. High RISC-recycling rates tend to reproduce the fully catalytic 
PTR profile, and low RISC- recycling rates approach the stoichiometric PTR 
profile. 

Recent experiments identified the average number of Argonautes in mam- 
malian cells to be of the order of 10^ proteins per cell^^. In human 293T cells, 
it is found that unregulated targets of 15000-25000 copies per cell are reduced 
approximately 5-fold by miR-223 reaching levels of 2700 copies per cell^. In 
order to test if our results and conclusions based on PTR in bacteria apply 
to other organisms, we have simulated PTR in mammals using the the above 
expression levels (data not shown). We find qualitatively similar results and 
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we conclude that our predictions hold also for mammalian PTR. 

In summary, our work elucidates how the molecular machinery behind 
PTR, and in particular competition of miRNAs and siRNAs during RISC- 
formation can influence the dynamics in target expression. We find two 
distinct regimes of PTR activity associated with low and high variation in 
the levels of target transcripts. It is possible for future work to construct a 
protocol in order to infer the mode of PTR, stoichiometric or catalytic, based 
on the variation of target expression levels, as well as elucidate the effects of 
Argonaute competition across different species of small RNA. 
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parameter 


definition 


value 


m 


small RNA number 




R 


target mRNA number 




P 


target protein number 




k 


small RNA production rate 


1/min 


dm 


small RNA degradation rate 


1/h 


kR 


mRNA transcription rate 


1/min 


dR 


mRNA basal degradation rate 


1/h 


dRm 


target dcstabilization rate 


1/h 


kp 


mRNA translation rate 


2/h 


dp 


protein degradation rate 


0.5/h 


kA 


RISC assembly rate 


(see Mo) 


dA 


RISC disassembly rate 


(see Mo) 


Ao 


Argonaute number 


500 


Mo 


d-A/kA 


1/600 - 600 


-Ro 


dA/dRm 


0.1 - 60 



Table 1: Definitions of model parameters. Typical values correspond to the activity 
of siRNAs in Salmonella ^^^lo^ Rates and concentrations are expressed in cell volume units. 
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Figure 1: PTR under high efHciency of RISC-formation. Fano factors (main 
graphs) and average molecular numbers per cell (insets) of the free small RNAs m (yel- 
low), loaded RISCs Am (blue) and free target transcript R (red) arc plotted as func- 
tions of the target gene transcription rate kfj under low competition PTR conditions 
{Mq — 1/600). Colored axis lines and axis labels indicate the corresponding quantity 
they refer to. Small RNAs have a fixed rate of production fc„i=l/min which determines 
the derepression threshold^. (A) Stoichiometric PTR. RISC and substrate degradation 
occurs in a short range around kn = km , and the transition from repressed to derepressed 
targets is sharp. (B) Catalytic PTR at 6-fold increased kfi range. Argonautes are titrated 
by the highly efficient substrate throughout the kfi range because RISCs are not degraded 
with the targets, which leads to strong target repression. 
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Figure 2: PTR under low efficiency of RISC-formation. Similar to Figure 1 but 
for low efficiency small RNAs (Mq = 600). (A) Stoichiometric PTR. Less efficient small 
RNAs are incorporated to lower numbers of Argonautes. The target repression to target 
derepression transition becomes smoother. (B) Catalytic PTR at 6-fold increased ku 
range. A distinctively different signature in target expression level variation is found. The 
strength in target repression is reduced compared to HE catalytic PTR, but the level of 
variation in target expression increases linearly with fc/j. 
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Figure 3: RISC-recycling determines the response of PTR to Argonaute 
scarcity conditions. The Fano factors of substrate small RNAs (yellow), loaded RISCs 
(blue), target transcripts (red) and target protein (gray) are plotted as functions of the 
Argonaute abundance Aq for LE small RNAs. (A) Stoichiometric PTR near the derepres- 
sion threshold fc/j — km — 1/min. More RISCs are formed as more Argonautes become 
available. However, RISC turnover increases also due to the joint RISC-target degradation, 
which results in an increase in the strength of target expression variation. (B) Catalytic 
PTR at fcfl = 12fc„i. Absence of joint RISC-target degradation renders the copy- number 
of small RNAs and the strength of variation of RISC copy- number independent of Aq. 
Lower numbers of available Argonautes lead to a reduced PTR-mediated target degrada- 
tion rate, which approaches the basal target degradation rate, resulting in increased levels 
of variation in target copy-numbers. 
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